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ABSTRACT 


The  objective  of  this  research  was  to  develop  beam  jitter  control  methods  for  a 
High  Energy  Laser  (HEL)  testbed.  The  first  step  was  to  characterize  the  new  HEL 
testbed  at  NFS.  This  included  determination  of  natural  frequencies  and  component 
models  which  were  used  to  create  a  Matlab/Simulink  model  of  the  testbed. 

Adaptive  filters  using  Filtered-X  Least  Mean  Squares  (FX-LMS)  and  Filtered-X 
Recursive  Least  Square  (FX-RLS)  were  then  implemented.  Disturbance  sources  included 
narrowband  mechanical  vibration  of  the  optical  bench  as  well  as  broadband  atmospheric 
turbulence  (simulation  only).  A  feedback  controller  with  adaptive  filter  (or  feedback- 
type  adaptive  filter)  was  applied  to  a  multi-rate  video  tracking  loop,  which  required 
precise  plant  identification  to  prevent  instabilities.  A  strapdown  type  Inertial  Reference 
Unit  (IRU)  system  was  investigated  using  a  reference  laser,  Position  Sensing  Device 
(PSD),  Fiber  Optic  Gyro  (FOG),  and  Fast  Steering  Mirror  (FSM).  The  controller  of  the 
strapdown  type  IRU  system  consisted  of  a  feedforward  gyro  signal  in  parallel  with  a 
feedback  PSD  signal.  A  Video  tracking  control  system  was  simulated  and  tested  on  the 
HEL  testbed.  The  IRU  strapdown  controller  was  simulated  and  preliminary  tests 
conducted.  Simulation  and  experimental  results  demonstrated  superior  performance  over 
classical  control  methods. 
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I. 


INTRODUCTION 


A.  MOTIVATION 

High  energy  laser  (HEL)  weapons  are  ready  for  some  of  the  most  challenging 
military  applications  in  future  battle  fields  since  speed  of  light  delivery  enables  the  war 
fighter  to  engage  very  distant  targets  immediately.  The  issues  of  the  technology  on  HEL 
systems  include  various  types  of  high  energy  laser  devices,  beam  control  systems, 
atmospheric  propagation,  and  target  lethality  issues.  Among  them,  precision  pointing  of 
laser  beam  and  high-bandwidth  rejection  of  jitters  produced  by  platform  vibrations  are 
the  key  technologies  in  the  emerging  fields  of  laser  communications  and  HEL  systems. 

Optical  beam  control  describes  the  centroid  shifting  of  a  laser  on  the  target,  and  is 
a  concern  of  engineers  and  scientists  working  with  lasers  and  electro  optical  systems. 
Platform  motion  and  optical  component  motion  causes  optical  jitter,  resulting  in  poor 
pointing  accuracy,  and  blurred  images.  Even  small  level  relative  motion  between  mirrors 
and  lenses  can  degrade  the  performance  of  precision  pointing  systems.  Sources 
contributing  to  optical  jitter  include  thermal  effects,  mechanical  vibration,  acoustics, 
static  and  dynamic  loading,  and  heating  and  cooling  systems. 

The  NPS  HEL  testbed  has  been  developed  to  support  research  environments  on 
the  precision  beam  control  technology  including  acquisition,  tracking,  and  pointing.  The 
testbed  incorporates  an  optical  table,  two  axis  gimbal,  high  speed  computers,  and  a 
variety  of  servo  components,  sensors,  optical  components,  and  software.  This  thesis 
describes  the  development  of  several  jitter  control  methods  and  their  implementation  on 
the  HEL  testbed. 

B.  THESIS  OVERVIEW 

Chapter  II  provides  background  information  relevant  to  this  thesis,  including 
information  on  HEL  systems,  beam  control  and  IRUs. 

Chapter  III  summarizes  the  HEL  testbed  here  at  NPS,  including  mathematical 
models  and  simulations. 
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Chapter  IV  provides  an  overview  of  the  various  control  algorithms  applied  to  the 
HEL  testbed. 

Chapter  V  describes  experimental  results. 
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II.  BACKGROUND  INFORMATION 


There  are  three  major  areas  that  provide  the  necessary  background  for  this  thesis, 
those  being  high  energy  laser  (HEL)  systems,  beam  control  methods — to  include  adaptive 
control  theory,  and  finally  inertial  reference  units  (IRU).  In  many  ways,  these  areas  are 
related  and  depend  upon  one  another  for  the  entire  system  to  work  properly.  Below  is  a 
brief  synopsis  of  these  areas. 

A.  HIGH  ENERGY  LASER  SYSTEMS 

Military  interest  in  HEL  systems  begin  shortly  after  the  invention  of  the  laser  in 
the  early  1960s.  It  is  easy  to  understand  such  interest,  what  military  would  not  desire  the 
ability  to  attack  an  enemy  at  light  speed,  incinerating  a  target  with  a  massive  burst  of 
photon  energy?  To  date,  three  major  demonstration  programs  in  the  field  of  directed 
energy  have  been  undertaken.  The  first  is  the  Army’s  Tactical  High  Energy  Laser 
(THEE).  This  was  a  megawatt  class  chemical  laser  using  deuterium  fluoride.  The  THEE 
program  proved  the  viability  of  HEL  systems  by  successfully  engaging  and  destroying 
several  Katyusha  rockets.  Engagement  ranges  were  on  the  order  of  several  kilometers 
[!]• 

Currently,  the  DOD  is  developing  two  airborne  HEL  systems,  the  Advanced 
Tactical  Laser  (ATL),  and  the  Airborne  Laser  (ABE).  Both  utilize  a  chemical  oxygen- 
iodine  laser  (COIL),  with  ABL  taking  the  open  cycle  approach  while  ATL  the  closed 
cycle.  Special  Operations  Command  (SOCOM)  is  the  sponsor  of  ATL  (see  Figure  1), 
which  is  mounted  in  the  belly  of  a  C-130  aircraft  and  intended  to  engage  moving  targets 
at  approximately  ten  kilometers.  Since  ATL  is  to  be  utilized  in  a  tactical  environment, 
there  is  a  need  for  a  closed-cycle  COIL  system,  which  captures  waste  chemicals  and  thus 
reduces  chance  of  detection  by  the  enemy.  As  of  October  2009,  ATL  has  successfully 
engaged  a  moving  target,  though  only  produced  limited  damage  [2]. 
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Figure  1 .  Advanced  tactical  laser  [From  2] 

For  strategic  applications,  the  Missile  Defense  Agency  (MDA)  is  developing 
ABL,  Figure  2.  This  system  is  the  largest  and  most  complex  to  attempted  as  of  yet.  ABL 
will  utilize  a  megawatt  class  COIL  integrated  on  a  Boeing  747-400F.  Operationally, 
ABL  is  intended  to  operate  above  40,000  ft  against  theater  ballistic  missiles  (TBM), 
which  will  be  engaged  during  the  boost  phase  [1].  In  addition  to  its  megawatt  class 
COIL,  ABL  also  has  two  kilowatt-class  solid  state  lasers  for  target  tracking/illumination 
and  atmospheric  measurements.  A  large  telescope  within  the  nose  mounted  turret  is  used 
by  the  beam  and  fire  control  systems  to  focus  COIL  beam  onto  a  pressurized  section  of 
the  target.  To  date,  ABL  has  successfully  acquired  and  tracked  a  boost  missile 
demonstrating  an  in  flight  targeting  sequence.  Additionally,  ABL  successfully  test  fired 
the  COIL  on  August  18,  2009  using  an  onboard  calorimeter  that  measured  beam 
characteristics.  MDA  hopes  to  successfully  engage  and  destroy  a  target  missile  by  the 
end  of  2009  [3]. 


Figure  2.  Airborne  laser  [From  3] 
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No  matter  the  application,  strategic  or  tactical;  airborne,  ground  based  or  sea 
based,  all  HEL  systems  require  extensive  beam  control.  This  is  due  primarily  to  the 
presence  of  self  induced  platform  jitter  and  atmospheric  disturbances.  Integral  to  proper 
beam  control  is  the  use  of  adaptive  control  methods.  Additionally,  beam  control  systems 
require  a  stabilized  line  of  sight  (LOS),  which  can  be  supplied  by  an  inertial  reference 
unit  (IRU).  Adaptive  filter  and  IRU  implementations  will  be  discussed  in  more  detail  in 
Sections  B  and  C  of  this  chapter. 


B.  BEAM  CONTROL  METHODS 

An  HEL  system  must  compensate  for  time  varying  disturbances,  which  requires 
the  use  of  adaptive  control  methods  to  reduce  beam  jitter  [4].  Like  classical  time 
invariant  control  systems,  adaptive  control  can  be  applied  using  feedforward  or  feedback. 
For  a  feedforward  system  (Figure  3),  a  coherent  reference  input  is  detected  prior  to  the 
control  actuators.  Such  a  reference  input  does  not  exist  for  feedback  system;  see  Figure 
4.  Both  feedforward  and  feedback  systems  rely  on  error  measurements  taken  at  the 
target.  In  the  case  of  an  HEL  system,  the  error  sensor  is  a  position  sensor  (PSD)  [5]. 
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Figure  3.  Feedforward  adaptive  filter  [From  6] 


5 


m 


In  Figures  3  and  4,  d(n)  represents  the  disturbance,  e(n)  the  error  as  measured  at 
the  target,  r(n)  the  reference  input  and  y(n)  the  output  of  the  digital  filter,  which  ideally  is 
identical  to  the  disturbance.  The  digital  filter  used  for  the  HEL  testbed  is  a  finite  impulse 
response  (FIR)  type  filter  (Figure  5). 


Figure  5.  Structure  of  digital  filter  [From  6] 

There  are  two  inputs  to  the  FIR  filter:  a  vector  of  filter  weights,  w(n),  and 
the  reference  signal,  r(n),  both  consisting  of  M  +1  stages.  A  time  delay  of  one  step  for 
each  of  the  M+1  stages  is  used  to  form  the  vector  of  time  delayed  inputs 
r(n)  =  [r(n),r(n-l),...,r(n-M  +  l),r(n-M)Y  .  Output  of  the  filter  is  a  control  signal 
y(n),  which  is  the  inner  product  of  the  weight  vector 
w(n)  =  [wq («),  Wj («),  W2 («),...,  (n)Y  and  reference  signal  vector  [7] : 
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j(«)  =  w^(n)r(«)  (1) 

The  error  signal  e(n)  is  the  difference  between  the  disturbance  d(n)  and  the  controller 
output  Equation  (2).  Unfortunately,  it  is  not  quite  as  simple  as  subtracting  from 
d(n);  one  must  also  take  into  account  the  dynamics  of  the  actuator  itself,  since  the  control 
signal  much  pass  through  a  fast  steering  mirror  (FSM)  before  its  effect  is  sensed  at  the 
target.  In  Equation  (2),  these  secondary  plant  dynamics  are  accounted  for  by  the  presence 
of  s(n  )\1~\. 

e(n)  =  d(n)-s{n)y(n)  (2) 

A  block  diagram  of  Equations  (1)  and  (2)  can  be  seen  in  Figure  6. 


Figure  6.  Block  diagram  of  transversal  filter  [After  8] 


At  each  time  step,  the  weight  vector  w(n)  must  be  updated,  which  takes  place  in 
the  Adaptive  Algorithm  block  of  Figures  3  and  4.  As  with  any  controller  method,  the 
objective  is  the  minimize  the  target  error.  An  easy  way  to  quantify  this  is  to  minimize  the 
mean  square  of  the  error  signal,  ^  =  E[e^(n)]  [5].  There  are  two  primary  methods  of 
updating  minimizing  4^ ,  least-mean-square  (LMS)  and  recursive-least-square  (RES). 
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1.  Filtered-X  Least  Mean  Squares 

Of  the  two  primary  methods,  least-mean-square  (LMS)  is  simpler  and 
computationally  less  expensive.  LMS  uses  a  gradient  approach  to  determine  the  steepest 
descent  for  minimizing  ^ .  In  [9],  it  is  shown  that  the  weight  vector  can  be  computed 
based  on  the  following: 

w(« -I- 1)  =  w(«) -t- //r(«)e(«)  (3) 

Steps  of  size  //  are  taken  as  the  fdter  converges  to  an  optimal  weighting  vector.  A  small 
step  size  generally  results  in  slower  convergence  to  a  value  closer  to  optimum,  where  as  a 
larger  step  size  may  produce  oscillations  about  the  optimum  weightings. 

As  stated  earlier,  controller  output  pass  through  the  secondary  plant  (in  our  case 
an  FSM),  prior  to  its  effects  being  measured  at  the  target  sensor.  If  unaccounted  for, 
phase  and  gain  variation  will  occur  between  the  error  and  reference  signals.  By  placing  a 

copy  of  the  secondary  plant  transfer  function,  iS(z) ,  in  the  reference  signal  path  to  the 
weight  updating  algorithm  of  Equation  (3),  such  variations  can  be  avoided  [7]. 

r{n)  =  s{n)  *  r(n)  (4) 

This  is  Filter-X  method  (FX-LMS),  which  prevents  an  unstable  controller.  In  simple 
terms,  the  reference  signal,  r(n),  is  filtered  prior  to  the  weight  updating  algorithm  rather 
than  prior  to  the  transversal  fdter  [7]. 

2.  Filtered-X  Recursive  Least  Meau  Squares 

Another  method  of  updating  the  weighting  vector  is  the  use  of  recursive-lease- 
squares  (RLS),  which  provides  faster  convergence  and  lower  steady  state  errors  than 
LMS.  The  same  FIR  fdter  is  utilized  for  RLS  as  LMS,  the  difference  being  the  algorithm 
used  for  weight  updating.  The  RLS  algorithm  contains  a  cost  function,  which  maintains  a 

memory  of  errors  according  to  a  forgetting  factor  of  0<A<1  (there  is  no  such  cost 
function  (<^)  memory  for  LMS,  which  only  uses  the  current  error). 
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Where  LMS  attempts  to  minimize  the  mean  square  of  the  current  error,  the  RLS 
algorithm  minimizes  a  summation  of  the  square  of  all  errors  multiplied  by  the  forgetting 
factor,  X  [7]: 


=  (5) 

i=l 

The  purpose  of  the  forgetting  factor  is  to  more  heavily  weight  recent  data  so  that 
nonstationary  disturbances  can  be  accounted  for.  Development  of  the  RLS  computations 
requires  a  significant  amount  of  linear  algebra,  which  is  explained  detail  by  Kuo  [7]  and 
Haykin  [9],  Equations  (6),  (7),  and  (8)  summarize  the  FX-RLS  algorithm  for  updating 
the  weighting  vector  w(n)  at  each  time  step. 


l  +  X-Y(n)Q(n-l)r(n) 


(6) 


w(n)  =  w(«  - 1)  +  (n)e{n) 


(V) 


Q(n)  =  X-'Q(n  - 1)  -  A-‘k(n)F  (n)Q(n  - 1) 


(8) 


As  with  the  FX-LMS  algorithm,  the  reference  signal  is  filtered  through  the  secondary 
plant  dynamics,  as  indicated  by  r .  Another  advantage  of  the  RLS  methods  is  that  the 
weighting  values  will  approach  those  of  the  optimal  Wiener  filter  weightings,  which  is 
not  the  case  for  LMS  [7]. 


3.  Reference  Signal 

The  FX-LMS  and  FX-RLS  algorithms  have  been  applied  using  a  single-channel 
reference  signal  in  the  research  of  [10],[1 1],  and  [12].  Additionally,  in  [8],  two  methods 
of  incorporating  multiple  reference  signals  are  implemented.  In  that  case,  the  reference 
signals  consisted  of  an  accelerometer  and  PSD.  In  this  thesis,  multiple  reference  signals 
are  also  utilized,  those  being  a  gyro  signal  and  again  a  PSD. 


C.  INERTIAL  REFERENCE  UNITS 

The  HEL  testbed  at  NPS  contains  an  imaging  sensor  with  a  frame  update  rate  of 
200  Hz  for  the  narrow  field  of  view  (NFOV).  Due  to  the  limited  frame  rate,  the  system  is 
unable  to  effectively  reject  high  frequency  jitter.  Additionally,  the  even  if  it  had  a  high 
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frame  update  rate,  the  NFOV  imaging  sensor  would  be  unable  to  detect  rigid  mode 
vibration.  For  these  primary  reasons,  an  inertially  stabilized  line  of  site  (LSO)  is 
required.  In  order  to  achieve  a  stabilized  LOS,  some  type  of  inertial  reference  unit  (IRU) 
is  needed.  Such  a  unit  will  control  high  bandwidth  FSMs,  which  cancel  the  effects  of 
optical  jitter.  There  are  two  primary  IRU  approaches;  stabilized  inertial  sensors  and 
strapdown  inertial  sensors.  Both  stabilization  methods  inject  a  reference  beam  into  the 
telescope,  which  follows  the  same  optical  path  as  incoming  imagery  as  well  as  the 
outgoing  high  energy  laser  [13].  Research  to  date  on  the  HEL  testbed  has  been  with  a 
strapdown  version  of  the  IRU. 


1.  Stabilized  Inertial  Sensor 


A  stabilized  inertial  sensor  system  is  one  in  which  a  stabilized  reference  laser  is 
mounted  to  a  gimbaled  platform.  See  Figure  7  for  a  simple  block  diagram  of  a  stabilized 
IRU.  A  closed  loop  control  system  is  used  to  drive  actuators  on  the  platform,  thus 
creating  an  inertially  stabilized  output.  The  reference  laser,  also  attached  to  the  gimbaled 
platform,  is  then  able  to  supply  an  inertially  stabilized  reference  beam.  Since  this  beam  is 
injected  into  the  same  optical  path  as  both  incoming  imagery  and  outgoing  laser  energy, 
one  is  able  to  determine  optical  path  jitter.  An  FSM  can  then  be  used  to  cancel  the  effects 
of  such  jitter.  Additionally,  track  sensors  are  used  to  command  the  platform  inertial 
sensors  and  therefore  compensate  for  target  motion.  A  major  advantage  of  a  stabilized 
inertial  system  is  that  an  exact  transfer  function  representing  dynamics  between  the 
reference  laser  and  the  corrective  element  is  not  required  [13]. 
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Figure  7.  Block  diagram  of  stabilized  IRU 
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2.  Strapdown  Inertial  Sensor 


The  strapdown  inertial  sensor  setup  is  the  simpler  of  the  two  methods.  In  this 
system,  inertial  sensors  are  rigidly  mounted  to  the  HEL  system  and  measure  base  motion, 
see  Figure  8  for  a  simple  block  diagram  of  such  a  system.  In  addition,  the  reference  laser 
is  rigidly  attached  to  the  inertial  sensors  and  thus  provides  alignment  of  the  optical  path 
to  the  sensor  frame.  A  feedforward  signal  from  the  strapdown  sensors  is  used  to 
command  an  FSM,  which  cancels  the  effects  of  jitter  on  the  optical  path.  As  with  the 
stabilized  system,  track  sensor  measurements  are  used  to  compensate  for  target  motion. 
While  strapdown  systems  are  simpler  and  cheaper  to  build,  they  are  prone  to  error  caused 
by  inaccurate  transfer  function  knowledge  [13]. 
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Figure  8. 


Block  diagram  of  strapdown  IRU 
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III.  NFS  HEL  TESTBED 


The  objective  of  the  HEL  testbed  is  to  provide  a  research  environment  for 
development  of  new  technologies  related  to  target  acquisition,  pointing,  tracking  and 
jitter  control  (see  Figure  9).  In  this  chapter,  a  brief  overview  of  the  testbed  is  provided. 
For  more  detailed  information  about  the  NFS  HEL  testbed,  refer  to  [6],  which  is  an  in 
depth  report  on  its  characteristics. 


i 

Figure  9.  HEL  testbed 

A.  SYSTEM  CONFIGURATION 

HEL  testbed  is  composed  of  three  major  components;  host  computer,  target 
computer,  and  beam  control  system  (Figure  10).  The  host  computer  manages  system 
operation  modes  and  all  sub-systems  through  user  interfaces.  Real-time  code  is  executed 
on  the  target  computer,  which  directly  controls  the  beam  control  system.  This  system  is 
composed  of  three  feedback  loops:  two  position  control  loops  and  one  rate  control  loop, 
as  shown  in  Figure  1 1.  In  addition  to  these  three  major  components,  a  moving  target  with 
an  illuminated  light  source  and  corresponding  position  sensing  device  (PSD)  provides 
evaluation  of  system  performance. 
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Figure  10.  HEL  system  configuration 


Figure  1 1 .  Control  block  diagram 


1.  Host  Computer 

An  MS  Windows-based  personal  computer,  on  which  all  software  is  developed, 
compiled  and  debugged,  acts  as  the  host  computer.  An  Ethernet  connection  is  used  to 
download  final  object  code  to  the  target  computer.  The  host  computer  also  provides  a 
means  for  user  input.  Motion  commands  are  generated  by  a  two-axis  joystick,  which  also 
includes  system  operation  control  switches.  Additionally,  menu  driven  interfaces  are  also 
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implemented  to  control  the  wide  field  of  view  (WFOV)  and  narrow  field  of  view  (NFOV) 
video  parameters.  Modes  include  video  mode,  track  mode,  gate  size,  and  video  tracking 
algorithm.  Figure  12  is  an  example  of  the  WFOV  user  interface,  a  similar  one  is 
provided  for  the  NFOV. 


Figure  12.  WFOV  user  interface 

2.  Target  Computer 

The  target  computer  is  a  Vxworks-supported  CompactPCI-based  system 
consisting  of  four  3U  size  boards:  PowerPC  CompactPCI  processor  board  (IMP2A),  10 
Pack  Carrier  Board  (ACP8630),  Multifunction  CompactPCI  Board  (ACPC730),  and 
Counter  Timer  Board  (ACP484).  Real-time  code  is  executed  by  the  PowerPC  board, 
which  also  controls  all  subsystems.  A  frame  grabber  PMC  card  is  also  mounted  on  the 
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board  to  control  and  communicate  with  the  WFOV  and  NFOV  eameras,  whieh  are 
conneeted  by  eamera  link.  This  board  also  communieates  with  the  host  computer  via 
Ethernet,  from  which  it  downloads  software  code  and  reeeives  control  commands. 
Angular  rate  data  is  received  from  the  gyros  by  a  synehronous  interfaee,  whieh  is 
connected  to  a  PMC  module  on  the  10  Paek  carrier  board.  Monitoring  of  analog  input 
signals  is  aecomplished  with  the  multifunetion  board.  See  Figure  13  for  an  external 
interfaee  diagram;  the  characteristies  of  the  four  boards  in  the  target  eomputer  are 
summarized  in  Table  1. 


Figure  13.  External  interfaees  of  target  computer 
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Table  1 .  Target  eomputer  characteristics 


Board  type 

Model 

Characteristic 

PowerPC  board 

IMP2A 

•  1.4GHz  PowerPC  7448,  3U  CompactPCI  SBC 

•  256  Mb  SDRAM,  1Mb  on-chip  cache,  128Mb  flash 

•  2  Ethernet,  2  serial  ports,  4  bits  GPIO 

•  PCI-x  capable  PMC  slot 

lO  Pack  Carrier  board 

ACPC8630 

•  Carrier  for  Industrial  I/O  Pack  Mezzanine  board 

Multifunction  board 

ACPC730 

•  16  bit  ADC  :  16  differential  or  32  single  ended, 

lOOKHz  conversion  rate  (lOuS  conversion  time) 

•  16  bit  DAC  :  differential  type,  80.8KHz  conversion 
rate  (12.375uS  conversion  time) 

•  32  bit  Counter/Timer  :  waveform  generation,  event 
counting,  watchdog  timing,  pulse  width  and  period 
measurement 

•  16  Digital  Input/Output  channels 

Counter  Timer  board 

ACPC484 

•  Six  32  bit  multifunction  counter/timer  :  position 
measurement,  pulse  width  modulation,  watchdog 
timer,  event  counter,  frequency  measurement 

•  16  digital  input/output  channels 

3.  Beam  Control  System 

Components  of  the  beam  control  system  include  WFOV  and  NFOV  track  loops, 
alignment  and  interface  optics,  reference  laser  and  surrogate  HEL.  A  schematic  of  the 
beam  control  system  is  shown  in  Figure  14,  while  a  breakdown  list  is  shown  in  Figure  15. 
To  create  an  optical  path  ffom/to  the  laser  source,  target  and  sensors,  a  number  of 
mirrors,  lenses,  and  beam  splitters  are  mounted  on  the  optical  table.  When  originally 
constructed,  the  HEL  testbed  also  included  an  analog  auto  alignment  control  loop,  which 
continuously  detected  and  corrected  misalignment  between  the  reference  laser  and  a 
position  sensing  devise  (PSD)  [6].  As  part  of  the  thesis,  this  analog  loop  was  replaced  by 
a  digital  controller  using  adaptive  control  methods.  The  new  auto  alignment  control  loop 
will  be  discussed  in  more  detail  in  Chapter  IV. 
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Figure  14.  Schematic  of  beam  control  system 
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Figure  15.  Breakdown  list  of  beam  control  system 


a.  WFOV  Track  Loop 

Two  feedback  control  loops,  inner  and  outer  loop,  comprise  the  WFOV 
track  loop  as  seen  in  Figure  16.  The  inner  loop  is  a  rate  control  loop  composed  of  a 
gimbal,  power  amplifier,  controller,  and  servo  components  (gyro,  motor,  and  encoder). 
Primary  purposes  of  the  inner  loop  are  to  maintain  LOS  to  a  target  in  inertial  space  with 
respect  to  external  disturbance,  and  track  input  rate  commands  from  the  WFOV  tracker. 
The  WFOV  camera  and  a  video  tracking  algorithm  comprise  the  outer  loop,  which  is  a 
position  control  loop.  This  loop  calculates  the  error  between  LOS  and  target  center, 
sending  an  error  signal  to  the  rate  command  of  the  inner  control  loop  [6]. 
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Figure  16.  WFOV  track  loop 
b.  NFOV  Track  Loop 

The  NFOV  track  loop  is  a  position  control  loop.  It  is  composed  of  the 
NFOV  camera,  fast  steering  mirror  (FSM),  and  a  video  tracking  algorithm,  see  Figure  17. 
As  with  the  WFOV  track  loop,  the  NFOV  track  loop  calculates  error  between  the  LOS 
and  target.  In  this  case,  though,  the  errors  are  much  smaller  than  can  be  compensated  for 
by  the  WFOV,  thus  minimizing  pointing  error  between  the  target  and  LOS. 
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Figure  17.  NFOV  track  loop 


B.  SYSTEM  CHARACTERIZATION 

A  number  of  experiments  were  performed  to  determine  characteristics  of  the  HEL 
testbed  with  the  results  used  for  system  modeling.  Experiments  included  resonant 
frequency  test,  rate  loop  servo  bandwidth  and  stabilization  test,  FSM  test,  and  NFOV 
bandwidth  test.  See  Figure  18  for  test  configuration.  An  external  terminal  board  on  the 
target  computer  provided  an  interface  for  signals  between  the  beam  control  system  and 
target  computer,  in  addition  to  providing  input/output  test  points.  Test  signal  generation, 
data  storage  and  observation  of  signals  were  conducted  using  a  dynamic  signal  analyzer, 
data  acquisition  system,  and  oscilloscope. 
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Figure  18.  Configuration  for  system  characterization  experiments 

1.  WFOV  Control  Loop 

a.  Resonance  Frequency 

Random  signals  were  applied  to  the  power  amplifier  and  output  signals 
were  received  at  the  gyro  and  encoder,  respectively.  Test  input  and  output  points  are 
shown  in  Figure  19.  Data  analysis  flow  for  resonance  frequency  determination  is  shown 
in  Figure  20.  Resonance  frequencies  were  calculated  using  power  spectral  density 
analysis.  Prior  to  performing  the  analysis,  preprocessing  was  required  to  remove  constant 
drift  in  addition  to  outliers.  A  periodogram  for  a  data  sequence,  [Vj,X2,---,x^] ,  was 
chosen  for  power  spectral  density  computation  as  given  by  Equation  (9).  To  suppress 
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spectral  noise,  an  average  of  20  power  spectral  density  computations  was  conducted. 
Resonance  testing  was  performed  over  a  frequency  range  of  0-200  Hz  with  an  input 


magnitude  of  1-2V,  as  summarized  in  Table  2. 
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Figure  20.  Analysis  flow  for  resonance  determination 
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Table  2.  WFOV  resonance  testing  scenarios 


Dir. 

Input 

(Vpeak) 

Output 

Data 

Time 

Length  (sec) 

Measured  Data  Analysis 

Freq.  Range 
(Hz) 

Average 

No. 

EL 

1.0 

Gyro 

8 

0-100 

20 

1.0 

Gyro 

4 

0-200 

20 

1.0 

Encoder 

8 

0-100 

20 

1.0 

Encoder 

4 

0-200 

20 

AZ 

1.0 

Gyro 

8 

0-100 

20 

1.0 

Gyro 

4 

0-200 

20 

2.0 

Encoder 

8 

0-100 

20 

2.0 

Encoder 

4 

0-200 

20 
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Over  the  wide  frequency  range  used  in  testing,  a  number  of  resonance 
frequencies  were  measured,  as  summarized  by  Table  3.  The  lowest  resonance  frequency 
of  the  elevation  axis  is  higher  than  that  of  the  azimuth  axis,  since  the  elevation  axis 
gimbal  mechanism  is  smaller,  more  rigid  and  simpler.  Additionally,  the  gyroscopes 
contain  many  more  resonance  frequencies  than  the  encoders  since  they  detect  angular 
velocity  versus  only  angular  position  for  the  encoder. 


Table  3.  Summary  of  WFOV  resonance  frequency  testing 


Direction 

Output  signal 

Resonance  frequencies  (Hz) 

AZ 

Gyro 

2.6,  8.3,  11.4,  15.9,  36.6,  69.5,  101,  134,  147.7 

Encoder 

2.6,  8.2,  11.4,  101 

EL 

Gyro 

3.8,  10,  14.9,  43.3,  61.4,  74.5,  101,  145,  195 

Encoder 

3.8,  11.6,61.4,  101 

b.  Rate  Loop  Servo  Bandwidth 

Rate  loop  bandwidth  was  measured  by  applying  a  sweep  sine  signal  to  rate 
command  and  observing  output  at  the  gyro.  Test  scheme  is  shown  in  Figure  21,  and  data 
analysis  flow  in  Figure  22.  As  with  resonance  testing,  preprocessing  was  conducted  for 
trend  removal  and  outlier  rejection.  Input  voltage  was  0.02V  for  elevation  axis,  and 
0.05V  for  azimuth  axis.  Data  was  measured  from  l-50Hz  for  both  axes. 


Figure  21 .  Rate  loop  servo  bandwidth  test  schematic 
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Figure  22.  Rate  loop  servo  bandwidth  data  analysis  flow 


The  ultimate  goal  of  this  testing  was  to  determine  the  transfer  function  of 
the  rate  loop  servos.  Several  intermediate  calculations  were  required  to  make  this 
determination,  the  first  of  which  is  the  cross-correlation  sequence  as  given  by 

RJm)  =  E[x„,^-y:]  =  E[x„-y„_:]  (10) 

where  x^  and  y^  are  jointly  stationary  random  processes,  and  £'[•]  is  the  expected  value 


operator.  In  reality,  a  correlation  estimate  given  by  Equation  (11)  had  to  be  used  since 
only  a  finite  segment  of  one  realization  of  the  infinite  length  random  process  was 
available. 


f  N-m-\ 

EX .  ■ V  m>0 

n+m  y  n 

Kxy{m)  =  <  /2=0 


Ry. 


m  <0 


(11) 


Additionally,  coherence  function  calculation  was  also  required.  Coherence  shows  the 
portion  of  the  output  spectrum  related  to  the  input  spectrum  as  given  by  Equations  (12)- 
(14) 


=FFT(RJ 

G„=FFnRJ 


(12) 

(13) 

(14) 


where  is  the  cross  spectrum  and  G^*  is  its  complex  conjugate,  G„  is  the  input  power 


spectrum  and  G^  the  output  power  spectrum. 
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Coherence  is  an  indication  of  statistical  validity  of  frequency  response 
measurements,  and  is  measured  on  a  scale  of  0. 0-1.0,  with  a  value  of  1.0  corresponding 
to  perfect  coherence.  Reasons  for  a  coherence  value  less  than  one  include  poor 
resolution,  nonlinearities,  extraneous  noise  and  uncorrelated  input  signals.  Additionally, 
it  is  independent  of  the  shape  of  the  frequency  response  function  since  it  has  been 
normalized.  Finally,  the  frequency  response,  also  called  the  “Transfer  function,”  is 
calculated  based  on  a  ratio  of  the  cross  spectrum  to  input  power  spectrum  as  given  by 

=  ^  (15) 

Gyro  output  for  the  azimuth  axis  is  shown  in  Figure  23,  while  elevation 
axis  gyro  output  is  shown  in  Figure  24.  Significant  steady  state  errors  were  observed  for 
both  axes.  Transfer  function  results  are  shown  in  Figures  25  and  26.  The  3dB  bandwidth 
is  6Hz  for  the  azimuth  axis,  and  7Hz  for  elevation.  Additionally,  resonance  frequencies 
in  the  azimuth  axis  of  8.2Hz,  1 1.4Hz,  and  15.9Hz  caused  significant  reduction  in  tracking 
performance. 


Figure  23.  Azimuth  axis  gyro  output 
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Figure  24.  Elevation  axis  gyro  output 


Figure  25.  Azimuth  axis  transfer  funetion 


Figure  26.  Elevation  axis  transfer  function 
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c.  Rate  Loop  Stabilization 

By  applying  a  disturbance  to  the  input  of  the  power  amplifier  and 
measuring  torque  error,  one  is  able  to  determine  the  torque  rejection  characteristics. 
Disturbance  input  was  a  0.5Vpk  sweep  sine  signal  over  the  frequency  range  of  1-lOOHz. 
Test  scheme  can  be  seen  in  Figure  27. 


Figure  27.  Rate  loop  stabilization  test  scheme 


Error  output  in  the  time  domain  is  seen  in  Figures  28  and  29,  while 
frequency  domain  output  is  shown  in  Figures  30  and  31.  Test  results  show  the  torque 
rejection  ration  is  low  over  the  tested  frequency  range.  Additionally,  stabilization  and 
servo  tracking  performance  were  adversely  affected  by  resonance  frequencies. 


OUMSIaMOlM 


Figure  28.  Azimuth  error  output 
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Figure  29.  Elevation  error  output 


Figure  30.  Azimuth  rate  stabilization  transfer  function 
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Figure  3 1 .  Elevation  rate  stabilization  transfer  function 

2.  NFOV  Control  Loop 

a.  Fast  Steering  Mirror 

FSMs  have  been  used  for  a  number  of  years  for  target  acquisition, 
scanning  and  beam  steering.  A  push/pull  configured  voice  coil  drives  a  two  axis  mirror. 
This  setup  is  similar  to  a  speaker  coil,  with  the  difference  being  the  FSM  is  configured 
with  a  moving  magnet  versus  a  moving  coil  as  is  the  case  for  a  speaker.  The  mirror  is 
flexurally  suspended  and  includes  an  optical  sensor  as  well  as  a  local  feedback  system. 
Local  position  feedback  is  the  inner  loop  of  the  NFOV  control  loop.  FSM  dynamic 
characteristics  were  determined  by  applying  a  sweep  sine  signal  to  the  local  position 
input  command  and  taking  output  from  the  position  sensor,  see  Figure  32.  Frequency 
range  for  the  test  was  1-lOOOHz  with  an  input  voltage  of  O.SVpk.  Tests  were  conducted 
on  both  azimuth  and  elevation  axis.  Results  from  frequency  response  tests  are  shown  in 
Figures  33  and  34.  As  one  would  expect,  the  transfer  function  in  nearly  identical 
between  the  two  axes,  and  both  have  a  -3dB  bandwidth  of  360Hz. 
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Figure  32.  Test  points  for  FSM  testing 


Figure  33.  NFOV  azimuth  FSM  transfer  function 


Figure  34.  NFOV  elevation  FSM  transfer  function 
b.  NFOV  Bandwidth 

Ideally,  one  would  be  able  to  use  a  test  scheme  as  shown  in  Figure  35  to 

determine  the  frequency  response  of  the  NFOV  video  tracker.  Unfortunately,  in  reality  it 
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is  not  so  simple  due  to  numerous  difficulties  generating  sweep  sine  target  motion,  and  a 
lack  of  target  position  values.  An  alternative  scheme  using  two  FSMs,  as  seen  in  Figure 
36,  was  actually  used  for  bandwidth  testing.  A  sweep  sine  signal  is  applied  to  the 
position  input  of  the  first  FSM  and  the  output  signal  is  measured  as  the  position  output  of 
the  second  FSM.  Test  signals  magnitudes  were  0.5V  and  l.OV  over  a  frequency  range  of 
O.l-lOOHz,  as  summarized  in  Table  4. 


Figure  35.  Ideal  test  setup  for  NFOV  track  loop  testing 
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Table  4.  Scenario  for  NFOV  track  loop  test 


Input  (Sweep  sine) 

Fine  Track  FSM 
Output  Pin 

No. 

Dir. 

Disturbance 
FSM  Input  Pin 

Sweep  sine 
Magnitude 

Frequency 

1 

X 

x+(4) 

0.5V(=1.31mil) 

0.1-100  Hz 

x-Pos  (pin  14) 

2 

(AZ) 

X-  (3) 

1.0V(=2.62mil) 

0.1-100  Hz 

GND(22,16) 

3 

Y 

(EL) 

y+(13) 

y-(12) 

0.5V(=1.31mil) 

0.1-100  Hz 

y-Pos  (23) 
GND(22,16) 

4 

1.0V(=2.62mil) 

0.1-100  Hz 

The  transfer  function  for  each  axis  was  nearly  identical  regardless  of  input 
magnitude.  In  both  axes,  the  -3dB  bandwidth  was  13Hz,  see  Figures  37  and  38  for  test 
results. 


Fraquaney  RawonaatAZ  Axil) 


Figure  37.  NFOV  video  tracker  transfer  function,  azimuth  axis 
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C.  MATHEMATICAL  MODEL  AND  SIMULATION 

A  mathematical  model  is  a  good  tool  for  estimating  system  performance  and 
evaluating  the  performance  improvements  associated  with  new  control  algorithms.  Two 
methods  were  used  to  determine  mathematical  models  for  the  various  HEL  testbed 
subsystems;  experimentation  and  component  specification.  Specifically,  data  sheets  were 
the  basis  for  models  of  the  servo  components  to  include  gyros  and  power  amplifiers. 
Based  on  the  component  models  and  experimental  results,  local  control  loop  and  whole 
integrated  models  were  built  using  Matlab  Simulink. 

1.  Component  Modeling 

The  KVH  Industries  model  DSP-3000  (now  referred  to  as  the  DSP-3400  by  KVH 
industries)  gyros  are  used  to  measure  angular  rate  of  rotation,  which  is  then  integrated  for 
accurate  angle  measurements.  The  DSP-3000  is  a  single-axis  fiber  optic  gyro  outputting 
a  digital  signal  over  a  synchronous  serial  interface  at  the  rate  of  1000/sec.  Bandwidth  is 
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an  important  parameter  for  dynamic  control  systems.  Since  the  -3dB  bandwidth  is 

greater  than  400Hz,  the  mathematical  model  of  Equation  16  is  used  for  the  gyro. 

2 


'GYRO 


-,w„=2=^;r*/,/  =  400,^  =  0.707 


(16) 


+  2^w^s  +  w/ 

In  the  HEL  testbed,  the  power  amplifier  utilizes  its  dual  phase  command  mode. 
The  phase  inputs  are  sinusoidal  and  120°  out  of  phase  from  each  other,  while  the  third 
phase  is  internally  generated  by  the  amplifier.  This  configuration  was  chosen  since  it 
provides  smooth  motion  while  minimizing  motor  torque  ripple  proving  maximum 
accuracy.  The  amplifier  mathematical  model  can  be  thought  of  as  an  EPF,  which  has  a 
-3dB  bandwidth  of  2KHz,  and  is  expressed  by: 


'  Power -AMP 


-,w„  =2*;r*/,/  =  2000,4^  =  0.707 


(17) 


s'^  + 

Optics  in  Motion  model  OIM-101  FSM  is  used  on  the  HEL  testbed.  These 
mirrors  consist  of  a  one-inch  glass  with  a  user-replaceable  mirror/sub-mount  where  the 
mirror  is  hard  mounted  to  the  mirror  gimbal.  A  built-in  high  precision  optical  sensor 
monitors  mirror  angles  and  provides  feedback  information  to  the  internal  controller.  A 
scale  factor  of  lOV  =  1.5  degrees  over  a  range  of  ±10V  is  used  by  the  local  position 
sensor.  As  previously  discussed,  frequency  response  testing  showed  the  -3dB  bandwidth 
of  the  closed  position  loop  is  360Hz,  meaning  that  the  FSMs  can  be  modeled  as  second 
order  low  pass  filter  (EPF)  given  by  Equation  18. 


^FSM  '  2 


-,w„=2*;rV,/  =  360, ^  =  0.707 


(18) 


-h  -h 

A  number  of  limit  switches  are  included  in  the  testbed  to  protect  various 
components.  See  Table  5  for  a  summary  of  these  limit  switches. 


Table  5.  Limit  function  summary  [From  6] 


Item 

Limit  Value 

Unit 

Voltage  Limit 

±10 

Vdc 

Current  Limit 

±10 

A 

Torque  Limit 

±23.9 

N-m 

Gyro  Limit 

±23.9 

Deg/sec 
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Finally,  a  disturbance  input  model  was  also  included  in  the  Matlab/Simulink 
system  model.  For  simplicity,  an  existing  ground  fighting  vehicle  model  that  expresses 
the  disturbance  as  a  power  spectral  density  function  was  utilized,  see  Equation  (19)19. 
The  RMS  disturbance  value  was  83mil/sec. 


PSD^ 


(mil/sGcf 

(l  +  4y2)3/2‘ 


,0.25<f<50Hz 


(19) 


2.  WFOV  Control  Loop 

Using  the  components  models  developed  in  the  previous  section,  the  WFOV 
control  loop,  consisting  of  an  inner  and  outer  loop,  was  constructed  as  shown  in  Figure 
39.  The  inner  loop  is  a  rate  control  loop,  which  provides  stabilization  with  respect  to 
external  disturbances  and  tracking  functions  based  on  gyro  feedback.  The  outer  loop  is 
an  angular  control  loop  that  automatically  maintains  LOS  to  the  target  center  based  on 
WFOV  camera  feedback. 


Figure  39.  WFOV  Simulink  model 
a.  Rate  Control  Loop 

The  step  response  and  corresponding  transfer  function  can  be  seen  in 
Figures  40  through  43.  Bandwidth  of  the  Matlab/Simulink  model  matches  that  of 
experimental  results.  The  frequency  response  shows  some  deviation  from  experiment 
results,  both  at  low  frequencies  and  around  the  mechanical  resonant  frequency.  This  was 
due  to  the  assumption  of  a  simple  linear  model  for  the  gimbals,  which  does  not  account 
for  resonance  or  nonlinear  effects. 
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Figure  40.  Azimuth  axis  step  response 


Figure  41.  Azimuth  axis  transfer  function 


Figure  42.  Elevation  axis  step  response 


36 


Figure  43.  Elevation  axis  transfer  function 
b.  Position  Control  Loop 

Outer  position  loop  simulation  results  are  shown  in  Figures  44  through  47, 
while  Table  6  contains  a  summary  of  model  parameters  and  simulations  results. 


Figure  44.  Azimuth  axis  step  response 
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Figure  45.  Azimuth  axis  transfer  function 


Figure  47.  Elevation  axis  transfer  function 
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Table  6.  WFOV  simulation  parameters  and  results 


Items 

AZ  axis 

EL  axis 

Gimbal  Inertia 

1000  in-oz 

260  in-oz 

Rate  loop 

Compensator 

94.4  ( P-controller ) 

28.8  ( P-controller ) 

-3dB  BW 

6Hz 

7Hz 

WFOV 

position 

loop 

Compensator 
(  Pl-controller) 

0.0063^  +  0.01982 

0.0069185  +  0.02173 

-3dB  BW 

5.4Hz 

5.9Hz 

3.  NFOV  Control  Loop 

The  model  for  the  NFOV  control  loop  is  shown  in  Figure  48.  The  FSM  is 
modeled  as  a  second-order  system,  as  was  verified  through  experimentation. 
Additionally,  the  NFOV  control  loop  includes  an  integrator  type  compensator  (on  each 
axis)  as  summarized  in  Table  7.  Simulation  results  for  one  axis  are  shown  in  Figures  49 
and  50,  since  each  axis  has  identical  characteristics.  Finally,  simulation  bandwidth 
matched  experimental  results. 


Test  Out 


Figure  48.  NFOV  Simulink  model 


Figure  49.  NFOV  step  response 
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Figure  50.  NFOV  frequency  response 


Table  7.  Summary  of  NFOV  control  loop  model  parameters 


Items 

AZ  axis 

EL  axis 

Compensator 
( I-controller ) 

1.2388 

1.2388 

-3dB  BW 

13Hz 

13Hz 

4.  Integrated  Control  Model 

The  WFOV  and  NFOV  models  were  combined  to  create  an  integrated  model  of 
the  HEL  testbed  (Figure  51).  The  WFOV  block,  containing  an  inner  rate  control  loop 
and  outer  position  control  loop,  is  located  in  the  lower  portion  of  the  integrated  model. 
The  upper  portion  of  the  model  contains  the  NFOV  track  loop.  Input  to  the  NFOV  loop 
is  the  track  error  generated  by  the  WFOV  track  loop. 
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Figure  5 1 .  Integrated  HEL  model 

Two  types  of  simulations  were  conducted  with  the  integrated  model.  In  the  first 
test,  target  tracking  performance  of  the  model  was  evaluated,  see  Figures  52  through  55. 
Target  motion  was  a  sinusoidal  with  2-degree  peak  at  a  frequency  of  0.5Hz.  At  steady 
state,  the  WFOV  has  a  peak  track  error  of  approximately  2.5mrads,  and  the  NFOV  has 
less  than  O.lmrads  of  peak  error. 


nraMNonon 


F igure  52.  T arget  motion 
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Figure  53.  Rate  of  WFOV  LOS 


Figure  54.  WFOV  track  error 
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Figure  55.  NFOV  track  error 


42 


In  the  second  test,  disturbance  rejection  characteristics  of  the  testbed  were 
investigated,  see  Figures  56  through  59.  The  disturbance  was  a  band  limited  random 
signal  with  83mil/sec  RMS,  as  described  previously  and  shown  in  Figure  56.  Peak  error 
for  the  WFOV  track  loop  increased  to  approximately  40  mrads,  while  NFOV  track  error 
increased  to  approximately  10  mrads,  see  Figures  58  and  59. 


Figure  56.  Disturbance  input 


Figure  57.  Rate  of  WFOV  LOS 
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Figure  58.  WFOV  track  error 


Figure  59.  NFOV  track  error 
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IV.  CONTROL  ALGORITHM  IMPLEMENTATION 


Two  control  algorithms  were  implemented  on  the  HEL  testbed.  The  first  was  an 
adaptive  filter  using  the  NFOV  video  tracker  only.  The  second  was  an  adaptive  filter 
using  a  strapdown  IRU,  which  consisted  of  a  PSD,  azimuth  and  elevation  gyros  and  a 
reference  laser.  The  NFOV  Simulink  model  described  in  Chapter  III  was  expanded  to 
include  multiple  disturbance  sources  (mechanical  vibration,  atmospheric  turbulence,  etc.), 
as  well  as  incorporating  all  three  HEL  testbed  optical  paths;  incoming  target  image, 
reference  laser,  and  surrogate  kill  laser;  see  Figure  60.  The  Controller  and  Plant 
subsystem  of  Figure  60  actually  contains  three  separate  controllers,  shown  in  Figure  61. 
The  first  is  the  video  tracker  feedback  loop,  the  second  a  gyro  adaptive  feedforward 
controller,  and  the  final  the  PSD  adaptive  feedback  controller.  Both  FX-LMS  and  FX- 
RLS  algorithms  were  implemented  for  all  three  controllers.  Most  of  the  component 
models  developed  in  Chapter  III  were  used  for  this  new  model. 


Figure  60.  NFOV  Simulink  model  incorporating  multiple  disturbances 
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F igure  6 1 .  NF O V  control  loops 

A.  ADAPTIVE  FILTER  USING  NFOV  VIDEO  TRACKER 

Adaptive  control  was  implemented  on  the  NFOV  using  only  the  video  tracker. 
This  implementation  was  of  the  feedback  type,  as  described  in  Section  II.B,  see  Figure  4. 
The  adaptive  controller  is  in  parallel  with  a  fifth-order  median  filter  (used  to  reject 
erroneous  noise  seen  in  experiments),  and  whose  output  is  prior  to  the  existing 
proportional-plus-integral  (PI)  controller,  see  Figure  62.  Additionally,  both  FX-LMS  and 
FX-RLS  adaptive  algorithms  are  implemented. 
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Figure  62.  Video  tracking  control  loop 

1.  Simulation  Results 

NFOV  video  tracker  simulations  were  conducted  with  three  disturbance  sources. 
The  first  two  were  from  mechanical  vibration,  one  at  20  Hz  and  another  at  50  Hz.  The 
third  disturbance  was  a  0.5  Hz  sinusoidal  in  addition  to  broadband  white  noise.  This  was 
to  simulate  the  effects  of  atmospheric  distortion  and  target  motion.  Only  the  PI  controller 
was  active  for  the  first  two  seconds  of  all  simulations.  A  comparison  of  FX-LMS  and 
FX-RLS  performance  is  shown  in  Figures  63  through  65.  The  frequency  analysis  shown 
in  Figures  63  and  64  was  conducted  on  the  steady  state  period  of  9-12  seconds.  As 
expected,  there  is  a  significant  reduction  in  target  error  once  the  adaptive  filter  is 
engaged.  Both  LMS  and  RLS  algorithms  reduced  RMS  target  error  to  less  than  one 
quarter  of  that  produced  by  the  PI  controller  alone.  While  the  FX-RLS  algorithm 
produces  lower  steady  state  error,  it  has  the  drawback  of  significant  transient  errors.  The 
FX-LMS  algorithm  is  much  slower  to  converge.  Further  indication  that  FX-LMS  is 
slower  to  converge  can  be  seen  in  Figure  65,  which  illustrates  the  convergence  of  FX- 
LMS  and  FX-RLS  weights.  After  ten  seconds,  FX-RLS  weights  have  stabilized,  while 
FX-LMS  weights  continue  to  change. 
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Figure  63.  Target  error  using  NFOV  FX-LMS  controller 


Figure  64.  Target  error  using  NFOV  FX-RLS  controller 
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Figure  65.  Video  tracker  adaptive  fdter  weight  comparison 

B.  ADAPTIVE  FILTER  USING  IRU 

Beam  control  contribution  from  the  IRU  consists  of  two  adaptive  filters  in 
parallel,  the  adaptive  feedforward  control  using  gyros,  and  the  adaptive  feedback  control 
using  the  PSD  (Figure  61).  The  IRU  is  exclusively  used  to  detect  and  correct  platform 
vibration,  it  will  not  detect  any  atmospheric  distortion.  While  the  NFOV  tracker  also  can 
detect  vibration  (sees  it  as  target  motion),  significant  errors  can  be  introduced  by  the 
video  tracking  algorithm.  Specifically,  the  video  tracker  uses  a  centroid  computation  to 
track  the  image.  Since  the  image  is  not  always  clear,  the  NFOV  video  track  loop  may  see 
errors  in  centroid  calculations  as  target  motion/platform  vibration.  The  IRU  is  less 
susceptible  to  this  problem  since  it  uses  a  laser  which  travels  much  shorter  distances  than 
an  image  that  is  seen  by  the  NFOV  video  tracker.  Also,  the  optic  path  the  reference  laser 
follows  is  not  subjected  to  any  atmospheric  distortion  which  aid  in  focusing  it  on  the 
PSD. 

The  gyro  adaptive  controller  was  implemented  using  FX-RLS  only,  while  the 
PSD  controller  was  implemented  using  both  FX-LMS  and  FX-RLS.  The  gyro 
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feedforward  eontroller,  Figure  66,  uses  the  gyro  signal  as  a  referenee  and  estimates  the 
error.  Speeifieally,  the  error  is  estimated  as  the  difference  in  the  gyro  signal  multiplied 
by  the  inverse  of  the  transfer  function  between  the  gyro  and  disturbance,  and  the  fdter 
output  multiplied  by  the  transfer  function  of  the  FSM,  see  Equation  20.  The  PSD 
feedback  controller  (Figure  67)  is  similar  to  the  video  tracking  adaptive  fdter,  since  both 
estimate  the  disturbance  based  on  the  error  and  fdter  output.  As  in  the  hybrid 
feedforward  and  feedback  controller  of  [8],  the  PSD  feedback  controller  error  is  the  sum 
of  the  gyro  feedforward  output  and  the  error  from  the  PSD  itself  The  filter  command 
calculated  by  both  IRU  controllers  is  summed  and  sent  to  the  auto-alignment  FSM. 

error  =  G-^^^g-Gf^J 
where : 

g  =  gyro  signal 

^gyro  -  transfer  function  between  disturbance  and  gyro 
f  =  fdtered  output 

^fsm  ^  FSM  open  loop  transfer  function 


Figure  66.  Gyro  feedforward  adaptive  controller 
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Figure  67.  PSD  feedback  adaptive  controller 

1.  Simulation  Results 

Simulations  were  conducted  using  the  same  0.5  Hz,  20  Hz,  and  50  Hz  sinusoidal 
signals  along  with  broadband  white  noise.  Simulation  results  of  the  PI  controller  only,  PI 
plus  PSD  adaptive  feedback  and  PI  plus  PSD  adaptive  feedback  along  with  gyro 
feedforward  can  be  seen  in  Figures  68  and  69.  For  both  plots,  the  IRU  adaptive  filters 
are  active  after  4  seconds.  The  frequency  analysis  is  from  the  steady  state  period  of  9-12 
seconds.  Using  the  adaptive  feedback  controller  combined  with  the  feedforward  adaptive 
controller  using  the  gyros  results  in  an  RMS  error  of  0.083  for  the  LMS  algorithm  and 
0.073  for  the  RLS  algorithm.  Additionally,  the  feedforward  controller  is  able  to  reduce 
the  20  Hz  disturbance,  which  the  PSD  feedback  filter  cannot  effect  (since  it  is  unable  to 
detect  rigid  body  motion). 
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Figure  68.  IRU  effectiveness  using  FX-LMS  Algorithm 


Figure  69.  IRU  effectiveness  using  FX-RLS  algorithm 
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V.  EXPERIMENTAL  RESULTS 


The  adaptive  eontrollers  from  Chapter  IV  were  implemented  on  the  HEL  testbed 
and  experiments  were  conducted  to  evaluate  their  performance.  In  all  experiments,  the 
disturbance  is  a  0.5  Hz  sinusoidal  signal  sent  to  the  disturbance  FSM.  This  disturbance 
affects  all  three  optical  paths.  Disturbance  magnitude  was  800  mV.  Other  natural 
disturbance  sources  include  the  power  amplifiers,  gyros  and  gimbals.  It  is  assumed 
atmospheric  disturbances  are  negligible. 

A.  STABILITY  ISSUES 

When  the  FX-LMS  adaptive  filter  was  implemented  in  the  NFOV  tracking 
algorithm  of  the  HEL  testbed,  stability  issues  arose.  After  approximately  10  seconds,  the 
system  would  become  unstable  and  lose  track  of  the  incoming  target  image.  It  was 
determined  that  the  stability  issues  was  caused  by  an  incorrect  model  identification  of  the 
Pl-Control  closed  loop  system.  The  more  accurate  plant  model  is  5th  order  (first  order 
controller,  second  order  FSM  and  second  order  track  sensor).  Figure  70  illustrates  the 
differences  in  frequency  response  for  the  two  models.  The  old  plant  model  has  a 
significant  phase  difference  above  10  Hz,  which  is  the  most  likely  cause  of  instability. 
Simulations  were  also  conducted  with  the  updated  plant,  see  Figure  71. 
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Figure  70.  Bode  plot  for  original  and  new  plant  models 


Figure  7 1 .  Simulation  comparison  of  old  and  new  plant  models 
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B. 


ADAPTIVE  FILTER  USING  NFOV  VIDEO  TRACKER 


Once  the  stability  issues  of  Section  A  were  identified  and  corrected,  FX-LMS  and 
FX-RLS  algorithms  were  implemented  on  the  NFOV  track  loop.  FX-LMS  experimental 
results  are  shown  in  Figures  72  and  73.  In  this  case,  the  NFOV  adaptive  filter  is  active 
after  approximately  35  seconds.  At  the  disturbance  frequency  of  0.5  Hz,  NFOV  track 
error  was  reduced  by  approximately  18  dB.  FX-RLS  experimental  results  are  shown  in 
Figures  74  and  75,  and  for  this  test  the  adaptive  filter  is  active  after  approximately  35 
seconds.  NFOV  error  was  reduced  by  approximately  19  dB  at  the  0.5  Hz  disturbance 
frequency  when  using  the  FX-RLS  algorithm. 


Figure  72.  NFOV  track  error,  FX-LMS  algorithm 
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Figure  74.  NFOV  track  error,  FX-RLS  algorithm 


56 


C.  FEEDFORWARD  ADAPTIVE  CONTROL  USING  GYROS 

Prior  to  implementing  a  strapdown  type  IRU,  a  demonstration  was  conducted 
using  the  elevation  gimbal  gyro  in  a  feedforward  signal.  To  accomplish  this,  the  transfer 
function  between  the  gyro  and  the  FSM  was  required.  This  transfer  function  was 
developed  using  two  tests,  the  first  to  determine  the  transfer  function  between  the  error  at 
the  PSD  and  the  gyro  output  and  the  second  a  transfer  function  between  the  FSM  output 
and  the  error  at  the  PSD.  Disturbance  motion  in  both  tests  was  detected  by  the  PSD. 
This  transfer  function  is  summarized  in  equations  20  through  22  where  is  the  error 

measured  at  the  PSD,  is  the  transfer  function  between  the  gyro  and  the  PSD,  g  is  the 
gyro  signal,  is  the  transfer  function  between  the  FSM  and  PSD  and  u  is  the  control 

signal  to  the  FSM.  Testing  in  Chapter  III  already  determined  the  FSM  could  be  modeled 
as  a  second  order  LPF  (Equation  18).  A  system  identification  was  conducted  to 
determine  G^ ,  which  is  a  first  order  transfer  function  (Equation  23). 

y,s,  =  G,g  (21) 

y psd  ~  ~GfSM^  (22) 
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(23) 


(24) 


Since  this  was  only  a  demonstration  of  the  feasibility  of  using  gyro  feedforward  control, 
the  disturbance  was  limited  to  5  Hz.  The  gains  for  both  and  were  optimized  for 


5  Hz. 


Experimental  results  for  the  gyro  feedforward  controller  can  be  seen  in  Figures  76 
and  77.  The  5  Hz  disturbanee  was  turned  on  just  prior  to  3  seeonds,  and  the  feedforward 
controller  was  turned  on  after  6  seeonds.  It  should  be  noted  there  was  a  Foresight  error, 
which  caused  the  target  error  to  be  non-zero  prior  to  application  of  the  disturbanee.  By 
using  the  gyro  feedforward  eontroller,  target  error  at  the  disturbanee  frequency  of  5  Hz 
was  reduced  by  approximately  8  dB.  Calibration  of  the  system  could  produce  further 
error  reduetions  but  was  not  attempted  since  these  experiments  were  meant  to  simply 
demonstrate  the  feasibility  of  use  the  gyro  signal  in  a  feedforward  eontroller. 


Figure  76.  Time  domain  gyro  feedforward  results 
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Figure  77.  Frequency  domain  gyro  feedforward  results 
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VI.  CONCLUSION 


A.  SUMMARY 

While  not  directly  related  to,  the  research  of  this  thesis  drew  on  much  of  the  work 
in  beam  jitter  control  accomplished  previously  at  NPS  by  Professor  Brij  Agrawal,  Dr. 
Joseph  Watkins,  Brett  Bateman,  Michael  Beerer  and  Dr.  Hyungjoon  Yoon  [8],  [10],  [12]. 
In  their  research,  feedforward  and  feedback  adaptive  fdters  using  multiple  reference 
sources  were  implemented  on  the  jitter  control  testbed  (JCT).  The  purpose  of  this 
research  was  to  implement  those  beam  jitter  control  techniques  onto  the  HEL  testbed. 

A  basic  system  characterization  of  the  HEL  testbed  was  required  prior  to 
implementation  of  any  adaptive  controllers  were  implemented.  This  characterization 
resulted  in  the  development  of  Matlab  Simulink  models,  which  proved  extremely  useful 
in  developing  and  testing  new  control  methods.  Adaptive  fdters  using  the  FX-LMS  and 
FX-RLS  algorithms  were  simulated  for  the  NFOV  video  tracker.  Additionally, 
simulations  were  conducted  using  an  IRU  for  a  hybrid  adaptive  controller.  This 
controller  contained  a  feedback  adaptive  filter  using  a  signal  from  the  PSD  and  a 
feedforward  adaptive  fdter  using  a  signal  from  the  gyros. 

Both  FX-LMS  and  FX-RLS  algorithms  were  implemented  on  the  NFOV  video 
track  loop  of  the  HEL  testbed.  Experiments  results  show  significant  reduction  in  track 
error,  especially  at  the  disturbance  frequency.  As  a  prelude  to  implementation  of  the  IRU 
hybrid  adaptive  controller,  a  simple  gyro  feedforward  controller  was  implemented  on  the 
elevation  axis  of  the  HEL  testbed. 

B.  FUTURE  WORK 

In  Chapter  V,  a  gyro  feedforward  controller  was  added  to  the  HEL  testbed.  The 
next  step  is  to  fully  implement  the  strapdown  IRU  and  associated  adaptive  controllers, 
which  were  simulated  in  Chapter  IV.  Currently  though,  the  gyro  feedforward  controller 
uses  the  gimbal  gyro,  which  is  located  at  the  base  of  the  telescope.  A  better  solution  is  to 
add  additional  gyros  at  the  tip  of  the  telescope,  which  are  co-located  with  the  reference 
laser. 
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Additional  work  would  be  to  develop  a  more  aecurate  transfer  funetion  between 
the  gyros  and  the  FSM.  This  could  be  accomplished  by  performing  tests  in  which  the 
loop  is  closed  and  error  measurements  are  taken  from  a  target  PSD.  A  more  accurate 
transfer  function  would  allow  the  gyro  feedforward  controller  to  perform  better  over  a 
wider  range  of  frequencies,  rather  than  being  tuned  for  a  single  frequency,  as  is  currently 
the  case. 

Finally,  this  research  concentrated  solely  on  implementation  of  adaptive  filter 
controllers  for  tracking  a  static  target.  Actual  HEL  systems,  such  as  ABE  and  ATE  must 
also  have  the  ability  to  track  a  moving  target.  The  WFOV  mode  of  the  HEL  testbed 
allows  it  to  track  moving  targets,  but  the  ability  of  the  NFOV  track  loop  to  maintain  lock 
on  a  moving  target  is  limited.  At  the  present  time,  not  work  has  been  conducted  on 
combining  the  WFOV  and  NFOV  control  loops  would  could  increase  the  ability  of  the 
NFOV  to  maintain  lock  of  a  moving  target. 
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APPENDIX:  EQUIPMENT  SPECIFICATIONS 


A.  ON-TRAK  OT-301  POSITION  SENSING  DEVICE 


Table  8.  ON-TRAK  OT-301  PSD  specifications 


Transimpedance  Gain  (V/A) 
Four  input  channels 

4x10^  to  4x10^  (Six  Ranges) 

Input  Current  Range 

1.5  //A  to  1.5  mA 

Output  Voltage: 

-position  X,Y 
-Sum 

+/-  lOV 

0-6  V 

Zero  Offset  (Offset  Null) 

+/-  1  V  Each  Axis 

Calibration  Adjust 

+/-  10%  Full  Scale  (+/-  IV) 

Detector  Bias 

0  V,  +/-  5  V 

Frequency  Response 

DC  to  15  kFfz  (range  dependent) 

G1  16  kHz 

G2  16  kHz 

G3  5  kHz 

G4  1.25  kHz 

G5  310  Hz 

G6  80  Hz 

Output  Connectors 

BNC 

Input  Connectors 

Sub  Min.  DB9  female  receptacle 

Power  Requirements 

12V,  500  mA  (AC  Adapter) 

Dimensions 

5.57  X  1.52  X  6.00  inches 
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B.  IMPERX  IPX  VGA  210-L  VIDEO  CAMERA 


Table  9.  IPX  VGA  2 1 0-L  video  camera 


CCD  sensor 

KAI-0340S/D 

Pixel  size 

7.4  //m 

Black  rows  -  top 

4 

Buffer  rows  -  top 

4 

Active  rows  -  (V) 

480 

Buffer  rows  -  bottom 

4 

Black  rows  -  bottom 

0 

Dummy  pixels  -  left 

12 

Black  columns  -  left 

24 

Buffer  columns  -  left 

4 

Active  pixels  -  (H) 

640 

Buffer  columns  -  right 

4 

Black  columns  -  right 

24 

Dummy  pixels  -  right 

12 

Frame  rate  -  single 

120  fps 

Frame  rate  -  duale 

210  fps 

C.  KVH  DSP-3000  FIBER  OPTIC  GYRO 

KVH  Industries  DSP-3000  series  fiber  optic  gyro,  part  no  02-1222-02  digital,  1 
kHz  synchronous  gyro. 
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Table  2:  Product  Specifications 


1  Attribute 

Rating  I 

Performance 

Digital 

Analog 

Maximum  Input  Rate 

±375“/sec 

j  ±1007sec  j 

Scale  Factor 

Linearity  (room  temp) 

1000  ppm,  1o  of  full  scale, 
for  ±375‘’/sec 

500  ppm,  la  of  full  scale 

500  ppm,  la  of  full  scale, 
for  ±1507sec 

Temperature  Sensitivity 

500  ppm,  la 

500  ppm,  la 

Error  (full  rate  &  temp) 

1 500  ppm,  1  a 

1 000  ppm,  1  a 

Bias 

Offset  (room  temp) 

±20°/hr 

- 1 

±1007hr 

Stability  (room  temp)* 

17hr.  la 

37hr,  la 

Temperature  Sensitivity 
(<1'’C/min) 

S^/hr,  la 

207hr,  la 

Bandwidth  (3  dB) 

>40  Hz, 

1 00  Hz  asynchronous 

100  Hz 

>400  Hz, 

1000  Hz  synchronous  or 
asynchronous 

Update  Rate 

1 00/sec, 

1 00  Hz  asynchronous 

2000/sec 

1000/sec, 

1000  Hz  synchronous  or 
asynchronous 

Angle  Random  Walk  (noise)* 

4°/hr/VHz 

0.06677Vhr  ( 

6“/hr/VHz 

0.17v'hr 

Initialization  Time  (valid  data) 

<5  sec 

*  Bias  Stability  and  Angle  Random  Walk  determined  by  Allan  variance  method. 


Figure  78.  KVH  DSP-3000  fiber  optic  gyro  specification,  part  1 
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Table  2:  Product  Specifications  (Continued) 


1  Attribute 

Rating 

Electrical 

Digital 

1  Analog 

Input  Voltage 

+5  VDC±10% 

Power  Consumption 

3  watts  maximum  (2  watts  typical) 

Output 

Type 

38,400  Baud  (RS-232), 

100  Hz  asynchronous 

±2  VDC  full  scale, 
differential 

115,200  Baud  (RS-232), 
1000  Hz  asynchronous 

±1  VDC  full  scale, 
single-ended 

3.072  MHz  serial, 

1000  Hz  synchronous 

Format  (selectable) 

Rate,  Incremental  Angle,  or 
Integrated  Angle 

Rate 

1  Physical 

Dimensions 

3.5“  X  2.3"  X  1.3" 

(88.9  mm  x  58.42  mm  x  33.02  mm) 

Weight 

0.6  lbs  (0.27  kg) 

Environ  mental 

Operating  Temperature 

-40°F  to  +167T  (-40°C  to  +75°C) 

Storage  Temperature 

-58‘’F  to  +185°F  (-50°C  to  +85°C) 

Shock  (functional) 

Functional  sawtooth  40  g,  6-1 0  msec 

Random  Vibration 

20  to  2000  Hz,  8  g  rms,  operational 

MTBF 

>55,000  hours,  ground  mobile 

In  all  cases  while  operating,  input  rate  must  not  exceed  ±500^/second.  If 
the  input  rate  exceeds  ±500° /second,  physical  damage  will  not  occur  but 
the  output  data  will  become  unreliable,  regardless  of  the  validity  BIT. 


Figure  79.  KVH  DSP-3000  fiber  optic  gyro  specification,  part  2 
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D.  OPTICS  IN  MOTION  OIMIOI  FSM 


Table  10.  OIMIOI  FSM  specifications 


Mirror  size 

1  in  X  0.25  in 

Wave-front  quality 

1/10  wave  p-v 

Useable  aperture 

0.94  in 

Scale  factor 

10  V  =  1.5  deg  (26.2  mrad) 

Range 

+/-  10  V 

Power 

+/-  15  V,  1.5  A 
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